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1.  INTRODUCTION 


This  report  covers  the  work  done  by  PhotoMetrics  on  the  GL/OPA  CO2  lidar 
system  including  optics  and  electronics  design  and  construction,  data  acquisition 
hardware  and  software,  testing  and  operational  procedures,  system  modeling,  field 
support  and  preliminary  data  reduction  for  the  SABLE  89  vAscension  Isl.)  and  GABLE 
(Azores)  field  measurement  programs. 

A  block  diagram  of  the  lidar  illustrating  its  operation  is  given  in  figure  1.1. 
The  transmitted  radiation  is  a  train  of  near  Fourier  limited  pulses  with  a  wavelength  of 
10.6  ^m  and  derived  from  an  injection  locked  CO2  TEA  laser  operating  on  the  P(20) 
transition.  The  characteristics  of  this  laser  are  listed  in  table  1.1.  Reception  is  coaxial 
with  transmission  and  via  a  25x,  eccentric  pupil,  Mersenne  telescope/beam  expander. 
As  shown  in  figure  1.1,  the  transmitted  and  received  radiation  are  separated  by  a 
polarizing  beam  splitter  and  quarter  wave  retarder  combination.  The  transmitted  light  is 
therefore  circularly  polarized  with  the  receiver  sensitive  to  circularly  polarized 
backscatter  of  the  opposite  sense.  The  receiver  employs  a  heterodyne  detection  scheme 
in  which  the  Doppler  shifted  light  backscatlered  by  aerosols  is  superimposed  at  the 
detector  surface  with  that  from  a  local  oscillator  (LO)  which  operates  at  a  fixed 
frequency  offset  (40  MHz)  from  the  transmit  laser.  The  resulting  intermediate 
frequency  (IF)  beat  in  the  detector  current  is  amplified  and  processed  according  to  the 
desired  mode  of  operation.  In  "intensity"  mode,  the  photocurrent  power  within  a 
bandpass  (B)  about  the  carrier  frequency  is  recorded  versus  time  of  flight.  The 
resulting  waveform  is  proportional  to  the  received  backscattered  power  versus  range. 
In  "complex  amplitude"  mode  the  IF  waveform  is  transformed  electronically  to  a  two 
component  baseband  signal  which  is  processed  to  yield  both  signal  amplitude  and 
frequency  versus  time.  From  these,  the  received  power  can  be  obtained  as  well  as  the 
radial  component  of  the  wind  velocity.  In  the  initial  operation  and  testing  of  this  lidar 
only  the  "intensity"  operational  mode  has  been  used. 

In  the  following  sections,  we  discuss  the  modeling,  sensitivity,  and  calibration 
of  the  system  first  and  then  proceed  to  the  details  of  the  lidar  layout,  alignment,  and 
electronics. 
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Table  1.1 

CO2  Laser  Operating  Parameters 


wavelength 

10.6  fim 

energy/pulse 

100  mJ 

pulse  duration 

1-2  fjLsec 

rep.  rate 

10-100  Hz 

spatial  mode 

TEMOO 

e'^  beam  width 

3.2  mm 

spectral  width 

<  500  kHz 

spectral  jitter 

=  1  MHz 

Figure  1 . 1  Illustration  of  lidar  heterodyne  operation. 
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2.  SYSTEM  MODELiNG  CALCULATIONS 


For  the  evaluation  of  system  performance  as  well  as  calibration,  an  analysis  of 

expected  signal  response  versus  range  as  well  as  its  dependence  on  the  wavefront 

properties  of  the  transmitted  beam  and  the  LO  is  needed.  A  complete  derivation  of  the 
"lidar  equation"  for  this  system  is  given  below.  We  have  found  such  an  ab  initio 
treatment  necessary  to  fully  appreciate  the  importance  of  all  the  physical  assumptions 
made  as  well  as  to  provide  a  unified  framework  for  modifications  should  ♦he  system 
performance  be  deemed  to  be  insufficiently  approximated  by  the  typically  assumed 
idealized  conditions,  e.g.,  diffraction  limited  optical  performance  and  negligible 
turbulence.  Although  the  treatment  given  below  is  generally  equivalent  to  that  in  the 
literature  [2.1],  v'e  adopt  a  formalism  and  notation  for  the  wavefront  analysis  which 

simplifies  the  resulting  expr-^'ssions  and  makes  more  transparent  their  physical 

interpretation. 

2.1  Lcfinitions  and  Units 

For  convenience,  we  adopt  units  for  the  electric  field  (E)  such  that  its  square 
(defined  to  be  the  field  intensity)  has  units  of  W/cm^.  The  spatial  dependence  of  the 
field  for  several  types  of  sources  is  presented  below. 

(i)  Gaussian  Beam: 

The  radial  (p)  dependence  of  ihe  field  for  a  Gaussian  beam  in  a  plane  perpendicular  to 
the  propagation  direction  is  given  by 

('’.1.1)  E(p)  =  A  (2/tw^)‘^^  exp{-p^/w^}  exp{-ikp^/2R} 

where  w  is  the  so-called  e'^  beam  radius,  R  is  the  wavefront  curvature,  and  k=2T/X. 
I  A  p  is  the  beam  power  in  Watts.  If  z  is  the  distance  from  the  beam  waist  to  the 
plane,  then  the  ra^hi  s  w=w^  and  curvature  are  given  by  the  standard  formulas 

[2.2]; 


(2.1. 2a)  w/  =  Wq^  [1 -l-(2z/kwo^)^l 
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(2.1.2b) 


=  z  [1  +(kwo^/2z)^]. 


In  these  e  xpressions,  v.q  is  the  mimmum  beam  radius. 

(ii)  Isotropic  Point  Source: 

The  fie’u  at  a  distance  r  from  an  isotropic  point  source  is  given  by; 

(2.1.3)  E(r)  -  Ae‘‘''/r(4T)’^2 

where  |  A  |  ^  is  the  total  power  emitted  by  the  source  fW). 

(lii)  Aerosol  Backscatter: 

As  long  as  the  receiver  subtends  a  small  enough  angle  at  the  particle  so  that  the 
differential  backscatter  coefficient  is  constant  over  the  entrance  aperture,  the  field  due 
to  the  backscattered  radiation  from  a  single  aerosol  particle  can  be  taken  as  that  due  to 
an  isotropic  point  source  with 

('2.1.4)  I  A  p  =  4t  iQ  3cr/5fli8o- 

Here,  Iq  is  the  incident  radiation  intensity  at  the  particle  (W/cm^)  and  da/dOigQ  is  the 
differential  backscatter  cross-section  (cm^/sr).  The  atmospheric  backscatter  coefficient 
((3)  is  defined  loosely  to  be  the  average  backscatter  cross  section  times  the  aerosol 
density.  More  rigorously 

(2.1.5)  /3  =  J  da  N(a)  da/anfa) 

where  N(a)  is  the  density  of  aerosols  with  radius  a.  If  the  treatment  is  to  include 
irregularly  shaped  particles,  the  integral  over  a  is  assumed  to  implicitly  include 
averaging  over  particle  shape  as  well  as  orientation.  Moreover,  dcr/dO,  and  therefore 
(3,  is  implicitly  dependent  on  the  polarization  of  the  transmitted  beam  as  well  as  the 
polarization  state  to  which  the  receiver  is  sensitive  (see  section  3.4). 
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2.2  State  Vector  Formalism  and  Dirac  N( 

The  paraxial  approximation  to  the  wave  equation  is  valid  for  the  conditions 
under  which  this  as  well  as  most  lidar  systems  operate.  Briefly,  only  rays  originating 
near  and  with  small  inclinations  to  the  transmit/receive  (i.e.,  optic)  tixis  need  be 
treated.  The  paraxial  approximi..ion  utilizes  this  to  reduce  the  second  order  wave 
equation  for  the  elecinc  field  to  one  which  is  first  order  in  the  longitudinal  coordinate 
(i.e.,  the  one  which  measures  distance  along  the  optic  axis— z).  D.  Stoler  has  shown 
[2.3]  that  the  resulting  equation  for  the  propagation  of  the  field  is  equivalent  to  the 
Schroedinger  equation  with  th^  propagation  coordinate  z  playing  the  role  of  time. 

It  is  natural  then  to  employ  the  standard  Dirac  bra-ket  notation  [2.4]  and 
represent  the  field  in  the  x-y  plane  at  z  by  a  state  vector  |E(z^>.  The  coordinate 
representation  of  this  state,  E(x,z)  as  used  in  section  2.1,  is  related  to  this  state  by: 

E(x,/)  =  <x  I  E(z)>  , 

where  |x>  is  the  coordinate  state  vector  which  is  an  eigenvector  of  the  two 
dimensional  traiisverse  coordinate 

operator  x  =  (x,y),  <xl  is  the  adjoint  stale  to  |x>,  and  <A|b>  represents  the 
inner  product  of  two  state  vectors.  Since  this  coordinate  state  representation  is 
complete  and  orthogonal,  the  following  identities  hold; 

<  X  j  X  '  >  =  5'^Hx-x  ’ ) 

1  =  J  dx  I  X  >  <  X  I  . 

where  is  the  Dirac  delta  function.  Using  these  relations,  the  inner  product  of  the 

states  representing  two  fields  can  be  s.iown  to  be  the  overlap  integral  of  the  fields: 

<  I- 1  I  E.  >  --  <  E.i  1  {  )■  dx  I  X  ><  X  I  }  1  E2  > 

=  Jdx<Eijx><x|E2> 

=  5  dx  Ei*(x)E2(x). 


and  with  the  unit  definitions  above  <e|E>  is  power  of  the  beam  wavefront 
represented  by  |  E  >  . 

Propagation  of  a  wavefront  along  the  z-axis  is  effected  by  a  unitary  operator, 
P(z'  ,z)  whic‘'  is  formally  equivalent  to  the  evolution  operator  in  quantum  mechanics; 

|E(z')>  =P(z',z)  |E(z)>. 

For  pro;  agation  in  free  space  or  in  an  isotropic  and  homogeneous  media,  the 
propagator  is  spatially  stationary  and  given  by 

(2.2.1)  Po(z)  -  e'*"^  exp{-izp2/2k} 

where  p,  the  "momentum"  operator,  has  a  coordinate  space  representation; 
p  =  -iV  =  -i(5/5x,9/3y). 

Propagation  by  this  operator  can  be  shown  to  be  equivalent  to  the  Fresnel  diffraction 
integral 


E(x,z)  =  (1/iXz)  J  dx'  e‘*“  exp{ik;(x-x '  )^/2z}  E(x',0). 

and  also  to  yield  the  standard  Gaussian  beam  formulas  discussed  in  section  2.1. 

The  representation  in  this  formalism  of  the  fields  discussed  in  the  previous 
section  is  straightforward.  The  transverse  dependence  of  a  Gaussian  TEMqo  beam  at 
its  waist  is  identical  ‘o  the  ground  state  wave  function  for  a  two  dimensional  (x-y) 
quantum  mechanical  harmonic  o,scillator  with  unit  mass  and  a  frequency  of  2/wo^.  We 
represent  this  field  by  |  (X)>  indicating  that  both  oscillate  s  are  in  their  ground  state. 
Higher  order  transverse  modes  of  a  stable  cavity— TEM^n,— arc  similarly  related  to  the 
excited  states  of  the  oscillator  - 1  nm>  .  It  was  in  fact  for  the  modal  analysis  of  stable 
cavities  that  this  quantum  mechanical  formalism  was  originally  adopted  for  the 
treatment  of  electric  field  wavefront  propagation. 

Under  paraxial  conditions,  the  field  at  (x,z)  due  to  a  point  source  at  (xo,z=0) 
can  be  written 
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Ep(x,z)  =  A  exp{ik(x-xo)^/2z}/z(4T)'^^. 


This  can  be  obtained  from  the  Fresnel  integral  equation  above  with 


Ep(x'  ,0)  =  iXA 

Since  <xlx'  >  s5^^\x-x'),  the  state  of  the  field  in  the  plane  of  the  source  can  be 
taken  as 

(2.2.2)  1f^(z  =  0)>  =  [iXA/(4x)’^2j 

With  a  redefinition  of  A  as  discussed  in  section  2.1,  this  also  is  the  state  of  the  field 
resulting  from  backscatter  from  a  single  aerosol  particle. 

2.3  Heterodyne  Signal  and  Noise  Formulas 

The  total  electric  field  at  the  detector  surface  (E)  is  the  superposition  of  the  LO 
field  (Elo)  ^"cl  that  due  to  the  received  backscatter  (Eg): 

1e>  =  exp{-iwst}  |Es(‘^>>  +exp{-iuLot} 

where  the  superscript  indicates  evaluation  at  the  detector  surface  and  wlq  (^s)  ’S  the 
angular  frequency  of  the  signal  (LO)  field.  These  frequencies  are  related  to  the 
transmitted  radiation  frequency  (cjj)  by 


“LO  -  wj  +  “c 


and 


—  a>j  +  OJQ. 


7 


In  these  expressions,  is  the  offset  (carrier)  frequency  between  the  LO  and 
transmitted  light  and  (jq  is  the  Doppler  shift  given  by 

wd  -4tv/X 

where  v  is  the  radial  velocity  of  the  scattering  particle  measured  positively  away  from 
the  receiver. 

The  photocurrent  induced  by  this  field  is 
i  =  (jjq/hi/)  <E  I  E> 

where  tj  is  the  photodiode  quantum  efficiency,  q  is  the  electronic  charge,  and  hi'  is  the 
photon  energy.  This  can  be  rewritten  as 

‘  =  ‘LO  +  ie  +  iiF 

where  Ilq  Ob)  current  due  to  the  LO  (backscattered)  Field  and  i|p  is  the 

current  at  the  intermediate  frequency  cj|F=a)c-a)D: 

‘LO  =  (’jq/hi')  <EloI  Elo>  =  (Who  Plo 

ie  =  (i?q/hj^)  <  Es  1  Es  >  =  (W^O  Ps 

ijF  =  (Who  [<Es|  Elo>  exp{ico,Ft}  +  c.c]. 

Plo  (Ps)  -s  the  LO  (received  backscatter)  power  incident  on  the  detector  (i.e., 
<  Epo  I  Elo>  )  c.c.  stands  for  the  complex  conjugate  of  the  previous  expression  in 
brackets. 

In  conventional  lidar  receivers,  there  is  no  LO  field  and  the  signal  current-ig- 
which  is  proportional  to  the  power  of  the  backscattered  radiation  is  the  measured 
quantity;  no  frequency  information  is  obtained  (unless  spectral  filtering  of  the  received 
backscatter  is  performed).  An  optical  heterodyne  receiver  measures  the  IF  current  so 
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that  both  amplitude  and  frequency  information  are  obtainable;  the  dc  signal  component, 
ig,  is  usually  negligible.  The  IF  signal  power  at  wjp  is  then 

(2.3.1)  p,F  =  2  Rl  (Wh.)2  I  <Elo^‘*>|Es<‘^>>  |2 

where  Rl  is  the  resistive  load  on  the  photodiode— e.g.,  the  input  impedance  of  the 
preamp— and  is  assumed  small  compared  to  the  diode  output  impedance.  It  is 
convenient  for  this  and  the  following  derivations  to  employ  normalized  field  states 
I  e>  defined  by 


e>  =  |E>/<E|E>‘^2 


so  that 


<  €  I  f  >  =  1. 

In  terms  of  this  representation,  the  formula  for  the  IF  signal  power  becomes  more 
explicitly  dependent  on  the  LO  and  received  signal  power; 

(2.3.2)  p,F  =  2  Rl  (7,q/h.)2  Plq  Ps  I  <  1  I 

Fiom  this  expression,  it  is  seen  that  the  IF  signal  power  is  proportional  to  LO 
power.  Therefore  the  IF  power  and  hence  the  signal  to  noise  ratio  (SNR)  increases 
with  increasing  LO  power  until  the  shot  noise  due  to  the  dc  LO  current  (Ilo)  dominates 
all  other  sources  of  noise  (e.g.,  detector  and/or  added  preamplifier  noise).  The  shot 
noise  power  in  any  bandpass  B  is  given  by  [2.5]: 

(2.3.3)  psN  =  2qBRL  Ilo  =  2qBRL  Plq. 

The  ratio  of  pip  to  psN  yields  the  SNR  for  heterodyne  operation  in  the  shot  noise 
limited  regime: 

(2.3.3)  SNR  =  Ps  I  >  P/hi/B. 
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~  ’?HET  V  Ps/hi'B 


where  r?HET  ‘s  the  so-called  heterodyne  efficiency: 


’7HET  -  I  <  f >  P- 
2.4  Derivation  of  the  Lidar  Equation 

We  have  implicitly  assumed  in  the  previous  section  that  the  detector  surface  is 
large  enough  to  subsume  the  LO  and/or  the  signal  field.  If  this  were  not  the  case,  then 
the  overlap  <  |  >  should  be  replaced  by  <  I  |  Es^*^^  >  where  IIj 

is  a  projection  operator  onto  just  that  portion  of  the  x-y  plane  occupied  by  the  detector: 

nd  =  i  d  dx  I  X  ><  X  1 

and  J  d  indicates  that  the  integral  is  restricted  to  the  detector  surface.  The  heterodyne 
efficiency  then  becomes 

’7HET  =  I  <  I  nd  I  >  I  ^ 

=  M  ddx  €lo(x)*€s(x)  I 

If  either  Elo  or  Es  is  negligible  outside  the  area  occupied  by  the  detector  the 
integration  can  be  performed  over  the  entire  plane  and  the  former  result  obtains.  We 
assume  this  to  be  the  case. 

Under  these  conditions,  it  is  not  necessary  or  even  convenient  to  evaluate  the 
overlap  integral  at  the  detector  surface.  For  example,  the  fields  at  the  detector  can  be 
obtained  by  propagating  those  exiting  the  focusing  lens  to  the  detector  surface,  so  that 


=Po(zid)  |Es<'>> 


|Elo^‘')>  =Po(z,d)  |Elo^‘>> 
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where  Po(z)  is  the  propagation  operator  as  discussed  above,  zy  is  the  distance  from  the 
lens  to  the  detector,  and  the  1  superscript  indicates  the  fields  at  the  lens.  Since  the 
propagation  operator  is  unitary,  Po^Po=l  and 


<  I  >  =  {  <  Elo<‘>  I  Po(Zld)^}  {Po(Zld)  I  Es^‘>  >  } 

=  <Elo<‘)|Es<*>>, 


the  overlap  can  be  evaluated  after  the  lens.  Since  the  transformation  of  a  wavefront  by 
a  lens  consists  of  just  a  multiplication  by  a  transversely  dependent  phase  factor,  it  too 
can  be  represented  by  a  unitary  operator  so  that  the  above  evaluation  may  be  performed 
before  the  lens.  In  fact  except  for  element  transmission  losses,  propagation  through  the 
entire  transmit/ receive  optical  system  can  be  represented  by  a  single  unitary  operator  as 
long  as  there  is  negligible  vignetting  of  the  received  radiation  and/or  the  virtual  LO 
beam.  By  virtual  beam,  we  mean  that  obtained  by  propagating  the  LO  beam  backward 
through  the  beam  splitter  and  into  the  receiver  leg  as  if  it  were  arriving  coaxially  with 
the  received  backscatter  radiation.  This  is  sometimes  referred  to  as  the  backward 
propagated  local  oscillator  or  BPLO  [2.1].  Thus  we  can  write 

I  «Lo'‘"  I  EsW  >  P  =  Tr  I  I  Es<"  >  P 

where  the  superscript  a  indicates  evaluation  of  the  field  wavefront  at  the  telescope 
entrance  aperture  and  Tr  is  the  receiver  transmission,  i.e.,  the  product  of  the 
transmissivity/reflectivity  of  all  of  the  optical  elements  of  the  receiver.  Note  that  we 
have  chosen  to  employ  the  full  signal  field  in  this  expression  while  for  convenience 
only  the  normalized  virtual  LO  field  is  treated  and,  for  its  evaluation,  transmission 
losses  are  to  be  ignored. 

This  procedure  can  be  carried  further  with  the  fields  being  propagated  out  to  the 
actual  scattering  volumes.  To  illustrate  this,  consider  the  signal  due  to  a  point  scatterer 
at  a  range  Zp  and  a  transverse  position  Xp.  Propagating  the  fields  out  to  the  plane 
containing  this  particle,  we  obtain  for  the  overlap; 
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I  <  I  Es"**  >  P  =  Tr  T^dp)  I  <  I  Es<f>  >  1 2 

-  TRTA(Zp)(ao/dO)  I  <  <Lo''"  I  Xp>  1 2, 

In  this  expression,  TA(Zp)  is  the  atmospheric  transmission  from  the  lidar  out  to  the 
range  Zp  and  the  expression  for  the  field  due  to  a  point  scatterer  from  section  2.2  has 
been  used.  The  incident  intensity  at  the  particle,  Iq,  is  given  by 

k  =  I  <Xp|&r(P)>  |2 

where  |  >  is  the  field  at  range  Zp  due  to  the  transmitted  lidar  pulse.  Using  the 

atmospheric  (T^i^)  and  lidar  transmitter  transmission  (Tf)  as  well  as  the  TEA  laser 
output  power  (Pj),  we  can  write  this  as 

Iq  =  Tj  Ta(Zp)  Pt  1  <Xp  I  I 

so  that  the  expression  for  the  overlap  is  explicitly  symmetrica]  in  the  transmitted  pulse 
and  the  virtual  LO  fields: 


1  <  I  >  P  =Pt  Tr  Tx  TP(Zp)  ida/dQ) 

'  I  <«LO^^  I  XpXXp  I  ex^P^>  I 

Using  this  in  the  formulas  from  the  previous  section,  we  obtain  for  the  IF  signal  power 
(referred  to  the  input  of  the  preamp)  expected  from  backscatter  off  a  single  aerosol 
particle  at  range  Zp  and  transverse  position  Xpi 

PlF(Xp,Zp;t)  =  2RL(jjq/hi/)2  Plo  Px(t-2zp/c)  TrTxTa^(Zp) 

(a<r/80)  X2-  I  <CLo^P)  I  XpXXp  I  6x(P)>  I  2. 

Note  that  the  temporal  dependence  is  due  solely  to  that  of  the  laser  pulse  power  and  the 
range  (z)  and  transverse  spatial  dependence-other  than  that  due  to  atmospheric 
extinction-arises  from  the  term: 
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f(x,z)  s  I  <  CLO^P)  I  X  ><  X  I  £j(P)  >  I  2 

=  I  <  eLO^“^  I  P(z.Za)^  I  X  ><  X  I  P(z,Za)  |  €7^“^  >  |  ^ 

=  I  €lo(x.z)«t(x.?-)  I  ^ 

where  P(z,Za)  is  the  propagator  from  the  telescope  aperture  to  the  range  z  and  eLo(*>^) 
(  ex(x,z)  )  is  the  normalized  virtual  LO  (laser)  field  at  this  range.  We  do  not  employ 
the  free  space  propagator  Po(z-zJ  as  atmospheric  turbulence  may  significantly  alter  the 
propagated  wavefronts.  In  principle  the  effects  of  turbulence  are  to  be  included  in  the 
detinition  of  the  propagator. 

For  convenience,  we  define  the  calibration  constant  A: 

(2.4.1)  A  -  2RL(i?q/hi/)2  Plq  TrT^  X^. 

Then  in  terms  of  this  the  IF  signal  due  to  a  single  aerosol  is: 

(2.4.2)  P[F(Xp,Zp;t)  =  A  PT(t-2zp/c)  TA^(Zp)  {daIdU)  f(Xp,Zp) 
and  that  due  to  a  collection  of  aerosols  at  {(XpZj)}  is 


p,F(t)  =  A  Ej  PT(t-2z,/c)  TA^(Zi)  (aa/an)i  f(x„Zi). 


By  passing  to  the  limit 

Ej  (dald(i)i  -*  )  dx  J  dz  0(z), 

we  obtain  for  the  expected  signal  due  to  atmospheric  aerosol  backscatter: 

(2.4.3)  PiF(t)  =  A  )  dz  P7(t-2z/c)  Ta^(z)  0(z)  F(z), 

where 
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(2.4.4) 


F(z)  s  J  dx  f(x,z). 

=  J  dx  1  eLo(x,z)eT(*.z)  1 

Equation  2.4.3  is  essentially  the  lidar  equation 

If  backscatter  from  a  calibration  target  at  a  range  zj  is  being  measured,  equation 
2.4.3  can  be  used  with 

/3  =  p*  5(z-zj), 

where  p*  is  the  target  material’s  backscatter  reflectivity,  i.e.  the  fraction  of  incident 
light  scattered  by  the  target  at  180°  per  steradian  and  under  the  prevailing  incidence 
and  polarization  conditions.  The  expected  signal  power  for  hard  target  calibration  is 
then 

(2.4.5)  pc(t)  =  A  p*  P-r(t-2zT/c)  Ta^(zt)  F(zt). 

For  calibration  purposes,  it  is  more  appropriate  to  use  the  time  integral  of  this  power  or 
the  energy  of  the  IF  pulse  due  to  return  from  the  target  referred  to  the  preamp  input: 

(2.4.6)  ec  =  A  p*  ej  TA^fzj)  F(zj). 

where  ej  is  the  laser  pulse  energy.  With  ec  known,  the  constant  A  can  be  eliminated 
from  the  lidar  equation  if:  p*  is  adequately  known;  the  atmospheric  transmission  to  the 
target  and  back  is  measured  or  accurately  modeled;  and  the  overlap  integral  F(z)  is 
known.  Although  only  an  explicit  range  dependence  for  and  F(z)  is  shown,  these 
quantities  are  implicitly  path  dependent  and  their  values  even  at  the  same  range  can 
vary  markedly  with  zenith  angle.  For  the  atmospheric  transmission,  this  variation  with 
zenith  angle  arises  mostly  from  the  fall  off  with  altitude  of  the  atmospheric  water 
vapor  content  which  is  the  dominant  source  of  extinction  at  10.6  /xm  at  low  altitudes. 
The  variation  of  the  overlap  with  zenith  angle  is  probably  less  drastic  but  also  more 
uncertain  as  it  arises  from  the  altitude  dependence  of  the  atmospheric  turbulence. 
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2.5  Modeling  of  the  Lidar  Overlap 


The  function  F(z)  as  defined  in  equation  2.4.4  is  the  overlap  integral  of  the 
normalized  transmitted  laser  intensity  and  the  virtual  LO  intensity  at  range  z.  The 
former  defines  the  irradiated  area  at  this  range  and  the  latter  the  transverse  dependence 
of  the  receiver’s  sensitivity.  F(z)  has  units  of  inverse  area  and  can  be  qualitatively 
interpreted  as  the  inverse  of  the  area  sampled  by  the  lidar  system  at  range  z.  No 
coherent  effects  are  evidenced  in  its  definition  (eq.  2.4.4)  because  we  have  referred  the 
calculation  to  the  plane  of  the  scatterers  and  there  is  no  correlation  in  the  phases  of  the 
backscatter  radiation  arising  from  different  particles.  Thus  F(z)  is  only  dependent  on 
the  intensities  of  the  fields  in  this  plane.  As  seen  in  the  previous  section,  calculation  of 
the  overlap  in  any  other  plane,  e.g.  the  telescope  entrance  aperture  or  the  detector 
plane,  results  in  apparent  coherent  effects,  i.e.,  field  dependence.  This  results  from  the 
partial  coherence  in  the  back^cattered  field  arising  from  propagation  [2.6]. 

The  evaluation  of  F(z)  requires  two  steps:  (i)  model  the  transmitted  and  LO 
wavefronts  at,  say,  the  telescope  entrance/exit  aperture  and  (ii)  model  the  atmospheric 
propagation  as  a  function  of  range.  The  former  should  include  any  vignetting  or 
aberration  in  the  transmitter/ receiver  optics  and  the  latter  the  effects  of  turbulence  on 
the  propagated  wavefronts.  In  this  preliminary  study,  we  will  neglect  these  effects 
effectively  adopting  the  free  space  propagator  Po(z)  as  defined  in  equation  2.2.1. 

The  ideal  case  which  most  closely  approximates  the  characteristics  of  the 
GL/OPA  CO2  lidar  is  that  in  which  the  transmitter  and  LO  are  Gaussian  TEMqo 
modes.  We  shall  assume  for  the  initial  study  that  these  wavefronts  are  well  matched 
and  collimated  at  the  exit  aperture  of  the  telescope.  In  this  case,  a  simple  calculation 
gives: 

(2.5.1)  F(z)  =  l/xw^^ 

where  w^  is  the  Gaussian  beam  width  at  range  z  as  given  in  equation  2.1.2a  and  the 
beam  waist  radius,  wq,  to  be  used  is  the  e'^  radius  at  the  telescope  aperture.  Using 
equation  2. 1 .2a,  this  can  be  rewritten  as: 

(2.5.2)  F(z)  =  (I/two^)  [1  +(z/zol2]-> 
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where  zq  is  sometimes  referred  to  as  the  Fresnel  or  diffraction  distance  and  is  defined 
by 

(2.5.3)  Zq  =  kwo^/2  =  twq^/X. 

It  is  at  this  distance,  zq,  from  the  waist  that  the  width  of  a  Gaussian  beam  has  increased 
by  a  factor  J  2.  It  is  at  this  distance  that  the  ideal  lidar’s  sensitivity  has  fallen  by  a 
factor  of  two  from  its  near  field  (z  <  <  zq)  value.  For  ranges  such  that  z  >  >  zq,  the  far 
field,  the  inverse  square  dependence  holds,  i.e., 

F(z)  =  k^wo^Mxz^  forz>>zo. 
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3.  SENSITIVITY  CALCULATION 


In  this  section,  we  perform  an  initial  sensitivity  calculation  for  the  lidar  system. 
As  seen  in  the  previous  section,  these  calculations  rely  on  ascertaining  a  number  of 
system  parameters  as  well  as  atmospheric  conditions  including:  atmospheric  and  optical 
transmissions,  turbulence  and  aberration  effects  on  the  overlap  function,  target  material 
reflectivity,  and  transmit  as  "veil  as  LO  beam  characteristics.  We  discuss  our  estimates 
of  these  parameters/effects  in  the  follo.ving  subsections.  In  the  last  part  of  this  section, 
v  e  compare  the  calculated  sensitivity  to  measurements  of  target  returns  obtained  during 
GABLE. 

3. 1  Beam  Sizes  and  Overlaps 

Given  the  TEA  laser  cavity  length  (L=  96  cm)  and  output  coupler  curvature 
(R=  10  m),  the  e  ^  radius  of  the  TEMqo  mode  at  its  waist  (the  internal  grating)  can  be 
found  from: 


wg  =  :X(R/L-1)‘'2j1/2 

=  3.2  mm. 

In  the  present  configuration,  the  distance  along  the  TEA  beam  from  the  grating  to  the 
telescope  secondary  is  about  3.8  m.  Using  equation  2,1.2a,  we  calculate  the  width  of 
the  TEA  beam,  at  the  secondary  to  be  5.1  mm.  This  gives  an  expanded  (25x)  beam 
width  (wq)  of  about  5.0"  and  a  Fresn  I  length  (see  equation  2.5.3)  of  4.8  km.  These 
results  are  sufficient  to  calculate  the  overlap  function  F(z)  versus  range  (from  equation 
2.5.2)  for  the  ideal  conditions  discussed  in  section  2.5. 

3.2  Optics  Transmission 

Table  3.2.1  and  3.2.2  list  our  estimates  for  the  transmittance/reflectance  of  the 
transmitter/receiver  optical  components.  The  reflectance  of  the  relay  mirrors  (Newport 
ER.2  coating)  has  been  estimated  from  the  manufacturers  published  typical 
performance  curve.  The  reflectance  of  the  SiO  coated  A1  scanning  and  telescope 
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mirrors  has  been  estimated  from  measurements  made  on  similar  mirrors  [3.1].  (Note 
that  a  lower  reflectance  for  the  scanner  mirrors  is  expected  due  to  the  45°  incident 
angle  of  the  beam.)  Typical  manufacturer  (II-VI,  Inc.)  specifications  are  assumed  for 
the  thin  film  jxilarizer  (TFP),  Rhomb,  lens  and  beam  splitter.  We  have  been  somewhat 
conservative  in  assuming  a  minimum  of  2%  insertion  loss  for  these  elements.  We  feel, 
however,  that  the  results  obtained  for  the  transmitter  and  receiver  transmissions  are 
high  due  to  expected  environmental  contamination  of  the  scanner  mirrors  (and  possibly 
the  telescope  primary). 


TABLE  3.2.1  Transmitter  Optics 


Item  Descr. 

ref /trans 

#  of  elements 

relay  mirror  (ER.2) 

0.98 

4 

tele,  mirror  (Al-SiO) 

0.97 

2 

scan,  mirror  (Al-SiO) 

0.95 

3 

TFP 

0.95 

1 

Fresnel  Rhomb 

0.98 

1 

Total  Transmission 

Throughput  (Tx)  =0.69 

TABLE  3.2.2  Receiver  Optics 


Item  Descr. 

ref /trans 

#  of  elements 

scan,  mirror  (Al-SiO) 

0.95 

3 

tele,  mirror  (Al-SiO) 

0.97 

2 

Fresnel  Rhomb 

0.98 

1 

TFP 

0.975 

1 

relay  mirror  (ER.2) 

0.98 

2 

beam  splitter 

0.93 

1 

lens 

0.98 

1 

window 

0.98 

1 

Total  Reception  Throughput  (Tr)  =  0.66 

18 


3.3  Atmospheric  Transmission 

For  calibration  as  well  as  evaluation  of  system  performance,  it  is  necessary  to 
correct  for  the  atmospheric  transmission  from  the  LIDAR  to  the  target.  As  the 
extinction  of  10.6  ^m  radiation  in  the  lower  atmosphere  (i.e.,  below  the  inversion 
layer)  and  under  typical  conditions  is  predominantly  due  to  absorption  by  the  water 
vapor  continuum,  the  transmission  is  highly  dependent  on  the  relative  h  '^idity  and  to 
a  lesser  degree  on  the  temperature.  Therefore  for  accurate  calibrated  backscatter 
measurements  it  is  necessary  to  characterize  the  atmospheric  conditions  and  model  or 
otherwise  measure  the  transmission. 

For  the  present  purpose  of  preliminarily  evaluating  performance,  we  take 
conditions  typical  of  those  encountered  on  the  GABLE  field  operation.  Table  3.3.1 
lists  the  results  of  PCTRAN  (based  on  LOWTRAN6  [3.2])  calculations  for  the 
atmospheric  trans  Mission  at  10.6  t-im  over  a  1  km  horizontal  path  as  well  over  the 
round  trip  distance  to  the  target  employed  on  GABLE.  We  adopted  the  standard 
midlatitude  summer  (ms)  atmosphere  with  the  Maritime  23  km  visibility  aerosol  model. 
The  default  values  were  taken  for  all  the  parameters  except  relative  humidity.  The 
range  of  RH  in  this  table  is  representative  of  the  low  altitude  radiosonde  measurements 
made  during  GABLE.  (Also  included  for  comparison  are  the  atmospheric 
transmissions  for  the  P(20)  transition  wavelength  of  CO2  from  a  listing  by  McClaichey 
and  Selby  [3.3]  for  tropical  (trop)  and  midlatitude  winter  atmospheres.)  We  indicate 
by  the  asterisk  the  value  for  T^i^^  which  will  be  used  in  sensitivity  calculation  below. 

TABLE  3.3.1  Atmospheric  Transmission  at  10.6  nvn 


Model 

Absorption  Coef 

T(lkm) 

T(6. 8km) 

PCTRAN  ms:  85%  RH 

0.34  /km 

0.71 

0.097 

PCTRAN  ms:  90%  RH 

0.38  " 

0.68 

0.073  * 

PCTRAN  ms:  95%  RH 

0.42  ” 

0.66 

0.059 

REF  3 . 3  trop 

0.61  " 

0.54 

0.016 

mw 

0.096  •' 

_ 1 

0.91 

0.52 
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Low  resolution  models  for  atmospheric  tran«mission  (such  as  LOWTRAN) 
should  be  sufficient  for  this  calculation  as  the  dominant  source  of  extinction  is  a 
continuum  and  not  a  band.  However,  we  feel  that  a  careful  evaluation  of  the  accuracy 
and  range  of  validity  for  the  available  atmospheric  models  should  be  made  in  the  future 
as  well  as  the  adoption  and  possible  testing  of  a  standard  model. 

3.4  Target  Reflectivity 

The  target  reflectivity,  p*  as  discussed  in  section  2.4,  is  by  definition  the 
fraction  of  incident  light  per  steradian  backscattered  toward  the  receiver.  It  is 
dependent  on  the  angle  of  incidence  to  the  average  target  normal  and  polarization  state 
of  the  beam  as  well  as  the  polarization  mode  to  which  the  receiver  is  sensitive.  For  the 
GL/OPA  CO2  lidar,  the  transmitted  beam  is  circularly  polarized  (CP)  and  the  receiver 
is  sensitive  to  CP  light  of  the  opposite  sense  .  (Note:  the  backscatter  of  CP  light  from 
Rayleigh  and  spherical  Mie  particles  as  well  as  specular  reflection  results  in  such  a 
reversal  in  the  light  helicity.) 

Due  to  its  portability  and  durability,  the  target  material  adopted  was  blue 
polystyrene  foam  [3.4].  Unfortunately,  the  target  material  characterization  that  JPL 
supplied  [3.5]  is  insufficient  without  further  assumptions  to  determine  p*  for  a  receiver 
sensitive  to  CP  light.  Basically,  JPL  measured  the  linearly  polarized  light  (lioth  s  and 
p)  backscattered  from  a  target  sample  with  CP  light  incident  at  45°  to  the  surface 
normal.  The  results  are  listed  in  table  3.4.1. 

TABLE  3.4.1  Measured  Target  Material  p* 


Incident 

Detected  Polarization 

Polarization 

s{LP) 

p(LP) 

right  (RCP) 

2.59  X  10"2 

4.17  X  10"^ 

left  (LCP) 

2.77  X  10“2 

4.45  X  10"^ 

average 

2.68  X  10~2 

4.31  X  10"2 
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It  is  not  clear  why  measurements  v/ere  made  with  both  RCP  and  LCP  incident 
light.  Only  an  unusual  and  highly  structured  material  surface  would  show  any 
difference.  The  6.5%  discrepancy  indicated  in  the  table  is  probably  a  better  estimate  of 
the  imbalance  in  the  JPL  quarter-wave  retarding  optics  than  of  any  physically  preferred 
sense.  Strengthening  this  argument  is  the  fact  that  the  ratio  of  s  to  p  reflectivities  in 
table  3.4.1  for  RCP  incident  light  differs  by  only  0.2%  from  that  for  LCP  incident 
light.  We  adopt  for  the  following  treatment  the  averages  of  the  measured  RCP  and 
LCP  reflectivities  and  refer  to  them  for  convenience  as  and  p^*. 

The  state  of  polarization  of  the  backscattered  light  is  unambiguously  given  by 
the  Hermitian  polarization  density  matrix  p  [3.6]  (not  to  be  confused  with  the 
reflectivity  p*): 


P  = 


P.S.S 

Psp 

.  Pp. 

Ppp  . 

1  Psp 

=  (Pps)* 

sensitive  to  some  pure  polarization  state  is  found  from 
(3.4.1)  P{^P)  =  ^Pp^p^ 


where  i/-  is  a  two  dimensional  row  vector  (and  its  adjoint-a  column  vector)  which 
represents  the  polarization  state.  For  example,  we  can  choose  to  work  in  a  basis  such 
that  the  representation.';  for  the  linear  s  and  p  polarization  filterings  are: 


4's  =  11,0] 

=  [0,1]. 

Then  p^^  and  Ppp  (as  the  notation  was  chosen  to  indicate)  are  the  fractions  of  light 
detected  through  s  and  p  orientated  LP  filters  and,  for  the  light  backscattered  from  the 
target,  are  given  by: 

Pv,s  ~  Ps  ^^Ps  ^  Pp  ^ 

Ppp  “  Pp  ^(p.s  ^Pp  )• 
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In  this  basis,  the  vector  representation  for  a  CP  filter  (with  sense  designated  simply  by 
±)  is: 


=  [l.±i]/(2)i'2 

and  the  fraction  of  light  detected  through  such  a  filter  can  be  obtained  from  equation 
3.4.1: 


=  0.5[1±  i(pps-Psp)] 

=  0.5[1  +  2  y  rj  sin<^] 

where  we  have  defined 

Pjp  —  y  r}  e"^ 

V  ~  (Pss  ^pp) 

The  parameter  y  is  related  to  the  degree  of  polarization  of  the  light.  If  7=1,  then  the 
determinant  of  the  polarization  density  matrix  is  zero  and  the  light  is  in  a  pure 
polarization  state  (generally  elliptical).  On  the  other  hand,  light  described  by  a  matrix 
with  7=0  is  as  "unpolarized"  as  it  can  be  consistent  with  Pss^Ppp-  The  angle  <l> 
describes  the  phase  of  the  residual  polarization,  e.g.,  0  =  ±-ir/2  for  light  which  is 
elliptically  polarized  with  +  sense  and  with  axes  as  defined  by  the  s  and  p  directions, 
and  4>=0  OT  TT  for  LP  light. 

It  is  clear  from  this  discussion  that  the  measurements  made  by  JPL— which  tell 
us  only  Pss  Ppp"^''^  insufficient  to  fully  characterize  the  polarization  state  of  the 
backscattered  radiation,  i.e.  there  is  no  way  to  obtain  7  and  <i>  without  making 
assumptions.  We  first  make  the  physically  reasonable  assumption  that  the  residual 
polarization  of  the  backscattered  light  is  elliptically  polarized  along  the  axes  determined 
by  the  target  orientation  and  with  a  sense  opposite  to  that  of  the  incident  CP  light.  Any 
other  possibility  would  require  the  target  surface  to  have  a  preferred  direction  (other 
than  its  normal)  and/or  sense  (i.e.  left-handed  vs  right-handed).  Given  this  assumption, 
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the  fraction  of  detected  light  for  a  receiver  sensitive  to  CP  light  of  the  opposite  sense  to 
the  transmitter  is 

P  =  0.5  [1+2  7  t}] 

and  for  the  effective  reflectivity  p* 

P*  =  P  (Ps"  +Pp*)- 

A  further  assumption  is  now  needed  for  the  value  of  7  where  7=0  for  complete 
depolarization  of  the  scattered  light  and  7=1  for  no  depolarizing  effects.  Here  no 
present  guide  is  available  and  we  take  7=0.5.  Under  this  assumption  we  obtain  for  the 
reflectivUy  appropriate  for  the  polarization  configuration  of  the  lidar:  p* =5.02x10'^. 

3.5  System  Performance 

For  the  preliminary  evaluation  of  system  performance  we  will  compare  the 
measured  IF  pulse  energy  due  to  backscatter  from  the  target  to  that  predicted  by 
equation  2.4.6  and  the  above  estimates  of  system  and  atmospheric  characteristics.  The 
latter,  using  equation  2.4.1,  is  given  by 

~  2RL(r;q/hj^)^  Plo  TrTj  X“  p*  ep  Tys^^(zj)  F(zj). 

=  4.7xlO'‘^J 

where  we  have  collected  the  assumed  quantities  below  in  table  3.5.1. 

For  a  comparison  to  measurement,  we  took  the  series  of  target  return 
measurements  AU210(X)1-AU21(X)12  taken  on  8/21/90  during  the  GABLE  field 
program.  Each  of  these  data  runs  consisted  of  averaging  the  results  of  100  laser  shots. 
Using  methods  to  be  discussed  later  in  the  section  on  electronics  and  calibration 
procedures,  we  obtained  an  average  return  IF  pulse  energy  of  4.3x10'^^  J  for  these 
measurements-about  10  times  lower  than  the  estimate  above.  Although  a  factor  of  3-4 
could  within  reason  be  accounted  for  by  our  overestimating  the  optics/atmospheric 
transmission,  detector  quantum  efficiency,  and  target  reflectivity,  we  have  no  doubt 
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that  the  major  portion  of  this  discrepancy  is  due  to  the  overestimate  of  the  overlap 
function  F(z)  which  resulted  from  our  assumption  of  idealized  conditions.  In 
particular,  what  must  be  accounted  for  and/or  corrected  is:  the  match  of  LO  and  TEA 
laser  beam  shapes;  aberration  and  possible  vignetting  of  the  transmitted  and  virtual  LO 
beam  by  the  telescope  and  scanner  mirrors  as  well  as  by  atmospheric  turbulence. 

TABLE  3.5.1  Target  Return  Sensitivity  Paiameters 


load  resistance  (Rl) 

50  0 

quantum  efficiency  (jj) 

0.5 

LO  power  (Plo) 

140  mW 

receiver  transmission  (Tr) 

0.66 

transmitter  transmissionCT-i-) 

0.69 

reflectivity  (p*) 

5.02x10-' 

laser  pulse  energy  (ej) 

75  mJ 

atmospheric  transmission  (T/^^(zt)) 

7.3x10-2 

overlap  (F(z-i')) 

13.1  m-2 
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4,  COo  LASER  AND  TRANSCEIVER  OPTICS 


Figure  4.1  illustrates  the  optical  layout  of  the  lidar.  The  labeled  components 


arc  idciitiiicu  below  in  lable  4.1. 


Tbe  lidar  iz  constructed  on  two  levels  with  the 


telescope  and  transmit/receive  (T/R  or  lower)  deck  suspended  from  the  laser  (or  upper) 
deck.  The  circles  in  the  figure  about  mirrors  MIO  and  M12  indicate  the  position  of  2" 
diameter  holes  through  the  optical  table  which  allow  transfer  of  the  TEA  laser  transmit 
pulses  and  the  local  oscillator  (LO)  beam  to  the  lower  deck.  The  position  of  these 
laser  ports  as  well  as  mounting  screw  holes  for  the  laser  are  shown  more  exactly  in 
figure  4.2.  Not  shown  in  the  diagram  is  a  large  turning  mirror  which  directs  the 
emergent  beam  vertically  into  a  hemispherical  scanner. 


Table  4.1 

Optical  Components 


Ml-2 

beam  steering  assemblies 

M3-20 

mirrors 

Bl-4 

beam  splitter/combiners 

IRl-4 

irises 

BX 

4x  beam  expander 

VA 

variable  attenuator 

FA 

fixed  attenuator 

CH 

chopper 

AOM 

acousto-optical  frequency  shifter 

Wl-2 

intra-cavity  Brewster  windows 

G 

TEA  laser  grating 

OC 

TEA  laser  output  coupler 

LI 

lens 

TFP 

thin  film  polarizer 

RH 

Fresnel  quarter  wave  rhomb 

DET 

cryogenic  HgCdTe  detectors 
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Figure  4. 1  Schematic  of  Transmitter/Receiver  Optics. 
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Figure  4.2  Optical  Table  Modifications. 


The  TEA  laser  itself  is  an  LSI  PRF-150  which  was  modified  by  LSI  for 
injection-locked  operation.  Additional  modifications  were  made  to  some  laser 
component  mounts,  the  injection  and  AFC  pickoff  scheme,  and  the  cw  beam  steering 
optics.  These  changes  were  made  to;  allow  for  installation  on  a  2’x7’  optical  table; 
reduce  the  transmission  of  vibrations  from  the  laser  head  fans  to  the  cavity  optical 
components;  allow  for  the  use  of  a  vertically  polarized  cw  laser;  and  to  eliminate  the 
necessity  for  a  beam  splitter  in  the  TEA  laser  beam. 


4. 1  Initial  Laser  Alignments 

4.1.1  CW  Laser 

Procedures  for  the  internal  alignment  and  operation  of  the  cw  laser  are  given  in 
the  MPB  manual.  We  discuss  below  procedures  for  initial  installation  and  alignment  of 
the  derived  beams. 

The  working  height  on  the  upper  deck  is  3.125"  and  on  the  lower  deck  4.875". 
The  exception  to  this  is  the  cw  beam  which  emerges  from  the  MPB  laser  vertically 
polarized  and  at  a  table  height  of  5.25".  This  beam  is  indicated  by  a  dashed  line  in 
figure  4.1.  The  portions  of  this  beam  picked  off  by  the  beamsplitters  B1  (the  seed  for 
injection  locking)  and  B2  (the  LO  for  the  AFC  detector)  are  displaced  to  the  upper  deck 
working  height  by  the  assemblies  Ml -2.  Due  to  the  configuration  of  these  assemblies, 
the  polarization  of  these  deflected  beams  is  horizontal.  There  is  no  beam  steering 
element  for  the  portion  of  the  cw  beam  which  traverses  the  beam  splitter  (the  signal 
LO).  Because  of  this,  it  is  important  that  during  installation  of  the  cw  laser  on  the 
optical  bench  the  output  beam  be  checked  to  ensure  that  it  passes  over  the  center  of  the 
laser  port  at  MIO.  When  this  is  done  the  mount  for  MIO  should  be  adjusted  vertically 
so  that  the  signal  LO  beam  strikes  the  center  of  MIO  and  is  deflected  down  through  the 
laser  port.  Further  alignment  of  this  beam  to  the  signal  detector  wil'  be  discussed  later 
in  the  section  on  the  T/R  optical  alignment. 

For  initial  alignment  of  the  seed  beam,  mirrors  Ml  and  M3  should  be  adjusted 
so  that  this  beam  passes  through  the  center  of  the  AOM  and  strikes  near  the  center  of 
M4.  Then  with  the  power  (and  cooling  water)  to  the  AOM  turned  on,  the  AOM  mount 
should  be  rotated  until  a  deflected  emerging  beam  is  seen.  This  indicates  that  insertion 
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is  at  the  correct  Bragg  angle  for  a  frequency  snift  of  40  MHz.  Fine  adjustments  of  the 
AOM  angle  will  maximize  the  power  in  the  shifted  component.  The  location  of  the 
beam  on  M4  should  be  checked  again.  Due  to  the  deflection  by  the  AOM,  it  may  no 
longer  be  near  the  center  of  the  mirror.  The  above  process  should  be  iterated  until  the 
shifted  seed  beam  can  be  deflected  by  M4  to  M5. 

For  the  AFC/LO,  mirrors  M2  and  M9  should  be  adjusted  to  direct  this  beam 
through  the  center  of  the  beamsplitter  B3  and  to  the  center  of  the  AFC  detector’s 
entrance  window.  This  should  be  done  with  the  AFC  detector  cooled  and  blocked. 
Then  with  the  chopper  (CH)  turned  on  and  the  detector  unblocked,  adjust  M9  to 
maximize  the  chopped  detector  current.  This  will  ensure  that  the  AFC/LO  is  centered 
on  the  detector.  (A  discussion  of  detector  and  biasing  operation  and  procedures  is 
given  in  a  later  section  on  the  receiver  electronics.) 

4.1.2  TEA  User 

Procedures  for  the  internal  alignment  and  operation  of  the  TEA  laser  are  given 
in  the  LSI  manual.  (Injection  locked  operation  is  not  contained  in  the  LSI  manual  and 
will  be  discussed  in  the  next  section.)  A  number  of  checks  and  alignments  with  the 
laser  firing  are  suggested  to  simplify  the  seed  and  AFC  pickoff  alignments  discussed  in 
the  next  section.  First  check  that  the  output  from  the  laser  is  exiting  through  the  center 
of  the  output  coupler  (OC).  If  the  beam  is  significantly  off  center,  correct  it  by 
"walking"  the  laser  axis  using  the  grating  and  output  coupler  tilt  controls.  Also  check 
that  he  (low  energy)  zero  order  reflection  off  the  grating  is  on  the  same  vertical  level 
as  the  center  of  mirror  M7  and  position  this  mirror  to  intercept  this  "pickoff  and  direct 
it  towards  M8. 

When  the  laser  operation  is  considered  satisfactory,  a  HeNe  laser  beam  is 
aligned  with  the  cavity.  To  do  this  remove  M20  and  Mil  (which  are  kinematically 
mounted  to  allow  for  precise  replacement)  and,  with  the  laser  firing,  center  the  iris  IR2 
on  the  TEA  laser  beam.  (We  have  also  found  it  useful  for  the  T/R  alignment  to  allow 
the  beam  to  propagate  several  feet  past  IR2— as  much  as  is  possible— and  strike/mark  a 
fixed  target  of  thermal  paper.)  With  the  TEA  laser  turned  off,  place  a  temporary 
mirror  in  line  with  the  TEA  laser  and  on  the  far  side  of  IR2  (e.g.,  near  the  table  edge) 
and  direct  the  HeNe  beam  at  this  mirror  by  rotating  M19.  Using  the  tilt  controls  for 
this  temporary  mirror,  M19,  and/or  the  HeNe  laser  mount,  walk  the  HeNe  beam  until 
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it  passes  through  the  center  of  the  iris  (IR2)  and  the  output  coupler.  A  number  of 
reflections  (probably  three)  off  the  output  coupler  should  now  be  observable.  This 
series  of  reflections  is  due  to  the  wedge  of  the  output  coupler.  The  middle  (or  second 
brightest)  reflection  is  due  to  reflection  off  the  inner  surface.  Continue  to  walk  the 
HeNe  beam  until  both  it  and  this  reflection  pass  through  the  center  of  IR2.  When  this 
is  done,  it  should  be  possible  to  trace  the  reflection  all  the  way  back  to  the  exit  aperture 
of  the  HeNe  laser.  The  HeNe  beam  is  now  lined  up  with  the  TEA  laser  cavity  axis. 

4.2  Injection  Lxxking:  Seed  and  AFC  Alignment 

4.2.1  Seed  Beam  Alignment 

The  seed  laser  is  to  be  injected  into  the  cavity  by  reflection  off  the  Brewster 
window  Wl.  V/ith  the  HeNe  laser  aligned  as  discussed  above,  a  dim  though  visible 
amount  of  HeNe  light  should  be  observable  exiting  the  side  window  on  the  mount  for 
Wl.  Position  M6  to  intercept  this  beam  and  then  adjust  the  tilt  of  M6  and  M4  so  that 
the  HeNe  beam  and  the  seed  beam  overlap  at  the  surface  of  M5.  Then  adjust  the  tilt  of 
M5  so  that  the  HeNe  strikes  the  center  of  the  exit  aperture  of  the  AOM.  The  seed 
beam  should  now  be  aligned  with  the  cavity  axis. 

4.2.2  AFC  Pickoff  Alignment 

With  the  HeNe  laser  aligned  as  above,  multiple  diffracted  beams  should  be 
observable  from  the  TEA  grating.  A  single  dim  order  (the  zeroth)  should  be  striking 
M7  if  this  mirror  was  positioned  as  discussed  in  section  4.1.2.  With  the  fixed 
attenuator  FA  removed  and  the  cw  laser  on,  adjust  M7  and  M8  so  that  this  beam 
overlaps  the  AFC/LO  at  the  surface  of  B3.  Then  adjust  the  tilt  of  B3  to  overlap  these 
beams  at  the  surface  of  M9.  This  last  step  is  difficult  due  to  the  low  brightness  of  the 
HeNe  light  reflected  off  B3  and  must  be  executed  with  all  room  lights  off.  After  this 
step  is  complete  replace  the  fixed  attenuator  and  insert  a  stop  in  front  of  IR2  to  prevent 
the  TEA  laser  from  striking  the  iris. 

4.2.3  Injection  Lxicked  Operation 

When  the  above  alignments  are  complete  and  with  the  AFC  detector  cooled, 
attach  the  detector  output  to  the  NERC  antplifier  and  the  output  of  the  amplifier  to  both 
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the  pulse  frequency  counter  (PFC)  and  an  oscilloscope.  Turn  on  the  amplifier,  PFC, 
and  the  EMT  driver  with  the  latter  in  closed  loop/unlocked  mode.  Turn  on  the  TEA 
laser.  With  the  AFC/LO  blocked,  a  detector  pulse  on  the  order  of  several  hundred 
millivolts  should  be  observed.  Fine  adjustment  of  the  alignment  of  the  TEA  laser 
pickoff  can  now  be  done  by  utilizing  M8  and  B3  while  observing  the  detected  pulse. 
When  the  AFC/LO  is  unblocked,  beating  between  the  LO  and  the  TEA  pulse  should  be 
observable  with  the  frequency  of  this  beat  varying  as  the  position  of  the  output  coupler 
is  manually  changed  with  the  EMT  driver.  As  the  output  coupler  is  slowly  scanned,  a 
fXDsition  will  be  reached  where  the  AFC  signal  waveform  dramatically  changes:  the 
beat  frequency  will  stabilize  at  40  MHz  and  the  laser  firing  time  will  advance  by  100- 
200  nsec  and  exhibit  reduced  jitter.  At  this  cavity  length,  the  frequency  of  the  TEMqo 
mode  of  the  TEA  laser  is  matched  to  that  of  the  injected  seed  radiation  and  the  output  is 
injection  locked.  The  PFC  counter  should  indicate  40  MHz  with  pulse-to-pulse  jitter 
less  than  about  1  MHz.  Switch  the  EMT  driver  to  open  loop/locked  mode.  This  will 
compensate  for  thermal  cavity  length  drifts  and  low  frequency  vibrations. 

4.3  Transmit/ Receive  Optical  Alignment 

4.3.1  Initial  HeNe  Alignment 

With  mirror  Mil  removed,  insert  M20  and  use  this  along  with  M19  and/or  the 
HeNe  laser  mount  to  line  the  HeNe  beam  with  the  center  of  IR2  and  the  spot  made  on 
the  fixed  target  in  section  4.1.2.  This  lines  the  HeNe  beam  up  with  the  TEA  laser 
beam.  Replace  Mil  and  position  it  so  that  the  HeNe  spot  is  centered  on  the  mirror. 
Now  use  Mil,  M12,  M13,  and  M14  to  pass  the  HeNe  beam  through  the  center  of  the 
TFP,  Fresnel  rhomb,  and  iris  IR4.  This  will  require  a  number  of  iterations  to  correctly 
position  the  TFP  and  rhomb.  When  this  is  finished  observe  where  the  HeNe  beam 
strikes  near  the  secondary.  Using  M12  and  M14,  walk  the  beam  until  it  passes  through 
the  center  of  IR4  and  strikes  near  the  center  of  the  secondary  mirror.  It  is  difficult  to 
judge  this  centering  on  the  secondary  so  that  for  fine  adjustments  the  position  of  the 
expanded  spot  on  the  primary  should  be  checked.  Now  center  the  iris  IR3  on  the  HeNe 
beam.  This  allows  the  two  irises  IR3-4  to  be  used  for  future  alignments  to  the 
telescope  optic  axis. 
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A  temporary  mirror  is  mounted  inside  the  telescof)e  side  access  panel  which  can 
be  swung  done  to  intercept  the  beam  several  inches  past  IR4.  Insert  this  mirror  and 
adjust  it  to  retroreflect  the  HeNe  beam  back  through  the  irises.  With  a  good  alignment 
this  retroreflected  beam  can  be  traced  all  the  way  back  to  the  HeNe  laser  exit  aperture. 
A  second  beam  reflected  by  the  TFP  towards  M15  should  be  observable.  Use  the 
position  and  tilts  of  M15  and  M16  to  direct  this  beam  through  the  center  of  the  beam 
splitter  B4  and  the  lens  LI. 

4.3.2  LO  Alignment 

With  the  signal  detector  removed  and  the  attenuator  VA  set  for  100% 
throughput,  use  mirrors  MIO,  Ml 7,  and  M18  to  deflect  the  LO  through  the  center 
B4.  By  placing  a  thermal  card  out  past  B4  towards  the  rear  of  M14,  both  the  Lo 
which  is  transmitted  through  B4  (95%)  and  that  portion  of  the  HeNe  beam  which  is 
reflected  by  B4  should  be  observed.  These  should  be  coincident  along  the  6-7"  path 
length  observable.  For  coarse  adjustment  use  M 17  and  M18  to  redirect  the  LO.  For 
fine  adjustments,  use  the  tilt  control  for  B4.  When  this  is  done,  check  the  alignment  by 
placing  a  thermal  card  at  the  detector  pxrsition.  The  focused  HeNe  and  LO  spots  should 
be  observed  to  be  coincident. 

To  align  the  detector,  turn  the  attenuator  down  to  a  setting  of  about  10,  cool  and 
replace  the  detector,  and  turn  the  chopper  on.  Then  adjust  the  detector  and  lens 
position  (including  focus)  to  maximize  the  chopped  detector  current.  (A  discussion  of 
detector  and  biasing  operation  and  procedures  is  given  in  a  later  section  on  the  receiver 
electronics.) 

At  this  stage  the  LO  is  lined  up  with  the  detector  but  the  receive  optics  can  only 
be  considered  to  be  coarsely  aligned  to  the  telescope  axis.  We  use  the  fraction  of  the 
seed  laser  which  is  transmitted  by  the  TEA  cavity  for  the  final  alignment.  To  do  this, 
remove  M20.  Now  the  weak  transmitted  seed  beam  traverses  the  transmit  optics  and  is 
retroreflected  to  the  detector  by  the  mirror  after  IR4.  To  observe  the  signal  produced 
by  the  seed,  connect  the  output  of  the  detector  to  the  signal  30  dB  preamplifier  and  the 
amplifier  output  to  an  oscilloscope.  A  40  MHz  signal  should  be  observed  due  to  the 
beating  between  the  seed  and  the  LO.  Moreover,  the  amplitude  of  this  signal  should 
vary  as  the  TEA  cavity  length  is  changed  by  manually  controlling  the  EMT  driver. 
The  TEA  laser  cavity  is  effectively  operating  as  an  etalon.  With  the  EMT  set  to 
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produce  the  largest  beat  signal,  mirrors  M15  and  M16  should  be  adjusted  to  maximize 
the  signal. 

When  this  stage  is  reached,  the  incident  LO  power  on  the  detector  should  be 
slowly  increased  by  raising  the  setting  of  the  attenuator  VA.  As  this  is  done,  the  beat 
signal  level  will  be  observed  to  slowly  increase  up  to  a  point  and  then  rapidly  decrease 
thereafter.  The  reasons  for  this  will  be  discussed  in  a  later  section  on  receiver 
electronics.  In  principle  the  attenuator  should  be  set  for  maximum  signal;  however,  to 
ensure  against  small  upward  cw  laser  power  drifts  it  is  preferable  to  choose  an  LO 
power  slightly  smaller  than  that  which  delivers  maximum  signal. 

4.3.3  Final  T/R  Alignment 

Due  to  the  difficulty  in  achieving  near  perfect  alignment  between  the  TEA  and 
the  HeNe  beams  (4.3.1),  we  have  found  it  necessary  to  perform  a  final  iteration  of 
some  of  the  above  procedures  .  ensure  that  the  transmit  and  receive  axes  are  aligned 
well.  First  with  the  mirror  past  IR4  removed  and  the  signal  detector  blocked,  the  TEA 
laser  is  turned  on  and  its  alignment  to  irises  1R3-4  observed.  In  addition  a  handheld 
detector  can  be  used  to  find  the  center  of  the  expanded  beam  at  the  output  end  of  the 
telescope.  If  the  result  of  the  observations  is  that  the  TEA  beam  is  aligned  with  the 
telescope  axis,  then  the  alignment  of  the  lidar  Is  finished.  If,  on  the  other  hand,  "mall 
displacements  are  noted,  then  mirrors  M14  and  M12  are  used  to  correct  these.  This 
adjustment  effects  the  pointing  of  the  seed  beam  which  was  used  for  the  receive  axis 
alignments.  Therefore  the  final  stages  of  4.3.2  must  be  repeated.  Notably,  with  the 
TEA  laser  turned  off  and  M20  replaced,  mirrors  M20  and  M19  along  with  the  HeNe 
laser  mount  are  adjusted  to  line  the  HeNe  beam  up  with  the  centers  of  IR3  and  IR4. 
The  temporary  mirror  is  then  inserted  on  the  far  side  of  IR4  and  adjusted  for 
retroflection.  M20  is  then  removed  and  the  beat  signal  is  maximized  as  above. 
Removal  of  the  retro  mirror  completes  the  alignment  procedure. 
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5.  RECEIVER  ELECTRONnCS 


5. 1  Detector  and  Bias  Circuitry 

The  detectors  employed  for  both  the  signal  reception  and  AFC  systems  are 
NERC  HgCdTe  liquid  N2  cooled  photovoltaic  diodes.  At  10.6  microns,  the  quantum 
efficiency  of  these  detectors  is  about  50%  with  response  on  the  order  of  4  A/W,  a 
signal  bandwidth  on  the  order  of  100  MHz  when  reversed  biased,  and  detector  noise 
v:ry  nearly  Johnson  noise  limited.  Figure  5.1.1  shows  the  typical  biasing  scheme  for 
these  detectors.  In  this  figure  standard  use  of  uf  per  case  letters  to  denote  the  dc  values 
of  current  and  voltage  is  employed.  Rp  and  Cp  are  the  diode  junction  shunt  resistance 
and  capacitance.  The  representation  of  the  detector  operation  by  idealized  components 
is  only  approximate.  In  practice  there  is  a  significant  degree  of  non-linearity  in  the 
diode  IV  curve  and  Rp  is  usually  taken  to  correspond  to  the  differential  resistance  at 
zero  bias.  Moreover,  the  junction  capacitance  decreases  and  the  detector  bandwidth 
increases  as  tiie  reverse  bias  voltage  increases.  This  results  from  *’  e  increase  in  the 
width  of  the  junction  depletion  layer  v/ith  reverse  bias. 

For  the  idealized  detector  illostrated  n  figure  5.1.1,  the  junction  voltage  is 
given  by 


where 

Vg  =  bias  supply  voltage, 

Rg  s  bias  .esiiiance, 

Ilo  —  dc  photocurrent  induced  by  the  LO  (=  response  x  Plq)- 

With  no  LO  radiation  incident,  the  junction  voltage  is  determined  by  the  current  (Ip) 
flowing  through  the  diode  shunt  resistance.  As  the  intensity  of  the  LO  is  incicased,  the 
junction  current  increases  which  reduces  the  amount  of  b’as  current  passing 
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Figure  5.1.1  Bias  Scheme  for  HgCdTe  Detectors. 


through  the  shunt  resistance  and  thereby  reduces  the  junction  voltage.  When 
Ilo=Vb/Rb,  the  junction  voltage  is  zero  and  there  is  no  reverse  bias.  This  explains 
the  behavior  discussed  in  section  4.3.2.  As  the  LO  field  increases,  the  IF  signal  first 
increases  as  expected  and  then  dramatically  falls  off.  This  decrease  is  due  to  the 
increase  in  junction  capacitance  (and  resultant  decrease  in  detector  bandwidth)  as  the 
reverse  bias  approaches  zero.  This  effect  provides  a  limit  in  addition  to  the  damage 
threshold  for  the  amount  of  LO  which  can  be  employed  since  the  manufacturer  suggests 
an  upper  limit  of  about  1  mA  for  the  bias  current.  In  the  next  section,  we  discuss  the 
implications  of  this  on  obtainable  signal  to  noise  levels. 

To  check  the  operation  of  the  detectors  under  different  bias  conditions  and  to 
measure  the  chopped  signals  used  to  align  the  AFC  and  signal  LO  to  the  detectors,  we 
constructed  a  bias  test  box.  The  schematic  for  this  circuit  is  shown  in  figure  5.1.2. 
This  box  allows  for  a  selection  of  a  wide  range  of  bias  ’•esistances.  As  discussed  in  the 
next  section,  typical  operation  from  a  9  V  supply  is  with  a  9  kft  bias  resistance  (i.e.,  a 
maximum  bias  current  less  than  about  1  mA).  Figure  5.1.3  shows  the  use  of  this 
circuit  to  monitor  the  chopped  LO  signals  employed  to  align  the  AFC  and  receiver  LO 
beams  (as  discussed  in  section  4.1.1). 

5.2  Petector/Preamp  Circuits 

For  the  analysis  of  detector/preamplifier  noise,  the  simple  model  illustrated  in 
figure  5  2. 1  is  used.  We  adopt  for  the  sake  of  example  the  specifications  of  detector  S- 
10.  From  the  NERC  test  report  these  are: 

output  impedance  (Rp)  300  0 

responsivity  4.8  A/W. 

The  amplifier  employed  is  an  Anzac  AM-113  with  the  following  specifications: 


gain 

30  dB 

noise  figure 

2dB 

input  impedance 

50  n. 
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Illustration  of  the  use  of  the  bias  test  circuit  to  monitor  chopped  CW  signals. 
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The  detector  noise  current  has  both  a  thermal  and  a  nonthermal— bias  voltage 
dependent— component.  The  rms  Johnson  noise  current  per  unit  bandwidth  for  a 
detector  temperature  (Tq)  of  77°  is 

io  =  4kTD/Rp  =  3.8  pAA/Hz. 

As  the  nonthermal  noise  current  in  similar  detectors  has  been  reported  to  be  less  than 
this,  we  will  take  the  above  as  an  estimate  of  the  total  noise  current.  The  quantities  i^^ 
and  Vy^  are  current  and  voltage  noise  source  models  for  the  preamplifier  added  noise. 
The  matched  source  noise  factor,  F,  is  related  to  these  by  [5.1]: 

(va)^/Rl  +  RL(iA)^  =  4kTo(F-l) 

where  Tq  is  the  laboratory  temperature  and  Rl  is  the  amplifier  input  impedance. 

For  the  AM- 113  (and  for  most  similar  IF  amplifiers),  the  relative  magnitudes  of 
va  and  iA  are  not  specified  as  only  matched  source  usage  is  envisioned.  We  consider 
the  three  cases  where:  (i)  the  added  voltage  noise  dominates  (va>  >RpiA);  (ii)  the 
added  current  noise  dominates  (iA>>VA/RA)i  and  (iii)  matched  conditions  apply. 
Under  these  conditions  the  effective  total  noise  power  delivered  to  the  load  can  be 
found  from  figure  5.2.1  to  be: 

case  i  Nl  =  4kRL  [RpTp  +  RlTo(F-1)]/(Rl+Rp)2 

(7.2x10-22  W/v^Hz), 

case  ii  Nl  =  4kRL  [RpTp  +  Rp2To(F-l)/RJ/(RL+Rp)2 

(7.6x10-21  WA/Hz), 

case  iii  Nl  =  k  [Tu  +  To(F-l)] 

(3.5x10-21  WA/Hz). 

In  the  following  discussion,  we  adopt  the  largest  of  these  estimates  (case  ii). 

To  obtain  the  largest  possible  signal  to  noise  ratio  the  detector  should  be 
operated  in  the  regime  where  the  shot  noise  due  to  the  LO  is  greater  than  or  equal  to 
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the  detector  and  added  preamplifier  noise.  The  mean  square  shot  noise  current  per  unit 
bandwidth  is  given  by 


*SH^  -  2  q  Ilo 


where  q  is  the  electron  charge  and  Ilo  as  above  is  the  LO  induced  dc  photocurrent. 
The  mean  square  detector  noise  current  which  would  deliver  Nl  to  the  preamplifier  is 

iN^  =  (Rp+Rl)2  Nl/RlRp^  =2.2x10-22  a2/Hz. 

For  the  shot  noise  to  dominate  the  detector/amplifier  noise:  ign  ^  in-  and  solving  for 
Ilo- 


Ilo  —  >N^/2q  =  0.68  mA. 


This  requires  the  LO  incident  power  to  be  greater  than  about  0. 14  mW. 

The  operating  condition  of  the  signal  detector  was  investigated  under  conditions 
afforded  by  the  manufacturer’s  suggested  maximum  bias  current  (1  mA).  With  a  bias 
supply  voltage  (Vp)  of  about  9  V  and  a  bias  resistance  of  about  9kn,  the  LO  power  was 
slowly  increased  until  a  maximum  was  observed  in  a  40  MHz  beat  signal.  For  this 
detector,  this  occurs  at  a  junction  voltage  (Vq)  of  about  280  mV.  The  bias  resistance 
current  (Ip)  at  this  setting  was  measured  (0.99  mA).  Then  the  local  oscillator  was 
blocked  and  the  bias  resistance  increased  until  the  same  junction  voltage  was  obtained. 
The  bias  current  under  these  conditions  was  measured  to  be  0.25  mA.  As  the  bias  and 
signal  current  are  related  by: 


Ib  =  Ilo  Ip(Vo), 


the  difference  between  these  currents  is  the  signal  current  induced  by  the  LO.  Thus 
under  these  bias  conditions  the  maximum  obtainable  signal  current  is  0.74  mA 
corresponding  to  an  incident  power  of  about  0.15  mW.  From  the  above  discussion  of 
expected  noise  levels,  these  results  indicate  that  detection  at  or  near  the  shot  noise  limit 
is  possible  under  these  conditions.  Indeed  qualitative  comparisons  of  the  noise  levels 
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observed  with  and  without  the  LO  indicated  an  increase  by  about  a  factor  two  in  noise 
power.  Note  that  for  such  a  comparison  to  be  valid,  the  noise  measurement  with  an 
incident  LO  beam  should  be  compared  to  the  noise  observed  without  incident  radiation 
and  the  bias  voltage  supply  or  resistance  altered  to  produce  the  same  junction  voltage. 
This  will  ensure  that  the  same  differential  source  impedance  and  capacitance  apply  in 
the  two  cases.  Future  evaluations  of  the  detection  system  should  include  quantitative 
measurements  of  the  obtainable  noise  levels  as  well  as  investigations  of  the  NERC 
detector  limitations  and  those  for  HgCdTe  detectors  available  from  other  sources.  A 
recent  study  [5.2]  indicates  a  lower  than  expected  heterodyne  quantum  efficiency  from 
similar  detectors.  This  may  be  related  to  bandwidth  considerations. 

Figure  5.2.2a  shows  the  schematic  for  the  preamplifier/bias  circuit  constructed 
for  the  signal  detection  system.  In  keeping  with  the  above  considerations,  it  will 
supply  about  1  mA  of  bias  current.  Figure  5.2.2b  shows  the  schematic  for  a  similar 
preamplifier/bias  circuit  constructed  earlier  and  along  lines  suggested  by  GE  [5.1]. 
Note  that  in  this  module,  the  bias  current  is  limited  to  a  value  less  than  about  0.5  mA 
and  according  to  the  above  discussion  will  not  be  capable  of  producing  shot  noise 
limited  operation.  This  module  can  be  easily  modified  to  deliver  a  higher  bias  current 
and  serve  as  a  spare. 

5.3  High  Gain  Amplifier  and  IF  Electronics 

Figures  5.3.1a-f  are  the  schematics  for  the  IF  amplifier  unit  which  we 
constructed  for  the  receiver.  This  unit  contains  two  separate  channels;  a  low  (<25dB) 
and  a  high  ( <  70  dB)  variable  gain  amplifier.  A  block  diagram  of  the  rf  module  layout 
is  shown  in  figure  5.3.1b.  A  complete  list  of  the  rf  modules  in  the  unit  (as  shown  in 
figure  5.3.1b)  is  given  in  table  5.3.1  with  manufacturer  specifications  in  Appendix  A. 
These  modules  are  all  connectorized  so  that  rerouting  and  substitution  is 
straightforward. 

We  have  included  in  this  unit  for  future  possible  use  the  log  amplifier  and 
frequency  discrimination  modules  which  were  delivered  to  GL/OPA  by  GE; 
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Schematic  of  Preamp/ Bias  circuit  (I B  =  1mA). 


44 


Figure  5.2.2b  Schematic  of  Preamp/Bias  circuit  (Ig  =  0.5mA). 


however,  these  have  not  been  used  or  tested  beyond  a  preliminary  check  of  their 
operability.  Moreover,  in  all  operations  to  date,  the  signal  split-offs  to  these  units  (the 
power  dividers  M12  and  M7)  have  been  bypassed  to  avoid  any  possible  interference. 
We  do  not  discuss  these  units  or  their  capability  in  this  report. 

Table  5.3.1 
RF  Module  List 


Module 

Make 

Model 

Description 

Ml 

MiniCircuits 

ZFSWA-2-46 

GaAs  SPDT  Switch 

M2 

EMF 

EMF-5006 

1 10  MHz  Oscillator 

M3 

ANZAC 

MD-140 

Mixer 

M4,12 

ANZAC 

THV-50 

Power  Divider 

M5 

K&L  ^Wave 

6LB3 1-70/5 

70  MHz  Filter,  5  MHz  BW 

6LB31-70/10 

70  MHz  Filter,  10  MHz  BW 

M6 

RHG 

LST60E13FP 

70  MHz  Log  Amplifier 

(13  MHz  BW) 

M7 

K&L  /xWave 

4LL32-90/I50 

Low  Pass  Filter,  90  MHz  cutoff 

M8 

RHG 

DT7020C 

70  MHz  Limiter/Discriminator 

M9 

RHG 

ICEVT7030C 

Linear  70  MHz  Amplifier 

(30  MHz  BW) 

Mil 

RHG 

ICF40E13FN 

Linear  40  MHz  Amplifier 

(13  MHz  BW) 

M13 

RHG 

DT40E13FN 

40  MHz  Limiter/Discriminator 

We  have  found  the  low  gain  amplifier  useful  for  utility  purposes  and  as  a  spare 
heterodyne  signal  preamplier  when  used  with  the  bias  test  box.  Its  13  MHz  bandpass 
(±6.5  MHz)  about  the  heterodyne  IF  frequency  (40  MHz)  may  make  this  amplifier 
useful  for  suppressing  out  of  band  transients  in  the  AFC  loop;  in  addition,  its  gain 
control  would  supply  the  capability  for  varying  the  AFC  signal  level  to  the  PFC.  The 
gain  of  this  amplifier  is  manually  adjustable  with  the  front  panel  knob  labeled  R9  in 
figure  5.3.1c. 
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5.3.  lb  Amplifier  unit;  RF  module  layout. 
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Figure  5.3.  Ic  Amplifier  unit:  front  panel  schematic. 
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Figure  5.3.  Id  Amplifier  unit;  rear  panel  layout. 


r-ii  vpc 
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Figure  5.3.  le  Amplifier  unit:  switching  and  MGC  board  schematic. 


Figure  5.3.  If  Amplifier  unit;  power  supply  circuitry. 


The  high  gain  channel  was  constructed  for  use  with  the  preamplified  heterodyne 
IF  signal.  This  unit  was  specified  to  cover  up  to  a  20  MHz  bandpass  about  the  IF 
center  frequency  (corresponding  to  a  Doppler  wind  velocity  range  of  ±50  m/s).  To 
obtain  this  coverage  in  a  high  gain  linear  amplifier  it  was  necessary  to  first  upconvert 
the  signal  center  frequency  to  70  MHz.  As  shown  in  figure  5.3.1b,  this  is  done  at  M3 
where  the  40  MHz  IF  signal  is  mixed  with  a  110  MHz  oscillator  signal.  The  resultant 
signal— centered  about  70  MHz-  is  fed  to  a  high  gain  amplifier  whose  bandpass  is  '’J 
MHz.  Because  the  detector,  added  preamplifier,  and  shot  noise  powers  are 
proportional  to  bandpass,  optional  filters  are  provided  to  reduce  the  amplifier  input 
bandpass  to  5  or  10  MHz.  These  are  especially  useful  for  increasing  the  signal  to  noise 
ratio  in  aerosol  backscatter  measurements  where  radial  wind  components  above  12.5  or 
25  m/s  are  not  expected. 

The  signal  source  for  this  high  gain  channel  can  be  switched  manually  or  in 
response  to  a  TTL  trigger  between  a  front  panel  (JIS:  source  A)  and  a  rear  panel  (J2: 
source  B)  connector.  With  the  front  panel  gate  mode  switch,  SW2,  in  the  manual 
position  (figure  5.3.1c),  the  source  switch,  SWl,  selects  between  A  and  B.  With  auto 
mode  and  A(B)  chosen,  source  A(B)  will  be  normally  selected  with  source  B(A)  gated. 
The  B(A)  selected  gate  period  begins  after  a  time  delay-adjustable  with  pot  R4  and 
range  switch  SW3-from  the  rising  edge  of  the  trigger  pulse.  The  duration  of  the  gate 
is  adjustable  using  pot  R6  and  range  switch  SW4.  This  switching  is  done  with  a  GaAs 
rf  switch  which  is  capable  of  switching  times  less  than  5  nsec.  In  practice,  we  have 
found  that  switching  transients  can  cause  the  downstream  high  gain  amplifier  to  go  into 
saturation  with  1-10  /xsec  recovery  times. 

We  have  also  built-in  a  triggerable  non-linear  gain  ramp  in  the  high  gain 
amplifier  channel.  With  switch  SW6  down  (ramp  off),  the  gain  is  manually  controlled 
with  knob  R7.  With  switch  SW6  up  and  SW7  down  (ramp  1  on,  ramp  2  off),  the 
initial  gain  is  set  with  R7  and  the  final  gain  with  knob  R8.  The  initial  gain  is  held  for  a 
time  after  the  rising  edge  of  the  trigger  pulse.  After  this  delay— adjustable  with  R12— 
the  gain  exponentially  "relaxes"  toward  its  final  value  with  a  time  constant  determined 
by  R14.  in  dual  ramp  mode— both  SW6  and  SW7  up— the  initial  gain  is  set  with  R7  and 
an  intermediate  gain  level  with  R8.  The  final  gain  is  the  maximum  attainable  with  the 
amplifier.  As  in  the  single  ramp  mode,  the  initial  gain  is  maintained  for  a  time  after 
the  trigger  pulse-controlled  with  R12-and  then  relaxes  to  the  intermediate  level  with  a 
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lime  constant  adjustable  with  R14.  After  a  second  time  delay  from  the  commencement 
of  this  first  ramp— adjustable  with  R 13— the  gain  ramp  time  constant  changes  to  one 
determined  by  R15  with  relaxation  towards  the  final— maximum— gain  level.  This 
operation  is  diagrammed  in  figure  5.3.2. 

Although  we  have  called  this  gain  ramping  exponential,  it  is  actually  the  gain 
control  voltage  which  is  varying  exponentially.  Since  the  gain  is  a  nonlinear  function 
of  this  voltage,  a  non-exponential  behavior  versus  time  of  the  gain  is  expected.  In 
addition,  rapid  gain  variations  are  limited  by  the  gain  control  response  time  of  the  RHG 
amplifier  to  about  a  10  fisec  time  constant.  For  setting  of  the  ramp  parameters  (as  well 
as  calibrating  gain  versus  time  or  range)  we  suggest  using  a  constant  power  input  to  the 
amplifier  at  40  MHz  and  monitoring  the  output  on  an  oscilloscope. 

5.4  Aerosol  Signal  Processing  and  Electronic  Calibration 

The  arrangement  used  for  aerosol  backscatter  measurements  during  the  GABLE 
field  program  is  illustrated  in  figure  5.4.1.  The  upconverted  (to  70  MHz)  and 
amplified  IF  signal  is  rectified  by  an  HP  rf  detector  diode.  The  resultant  signal  is 
digitized  and  averaged  over  a  selectable  number  of  lidar  shots  by  a  DSP  2112FS 
CAMAC  module  signal  averager.  The  averaged  waveform  is  passed  to  and  recorded 
by  the  AST386  computer  which  controls  the  CAMAC  era  ..  (On  SABLE  89,  a  similar 
configuration  was  used  but  with  the  Analogic  D6000  employed  to  digitize,  average, 
and  record  the  signal  waveforms.) 

The  rf  diode  is  roughly  a  square  law  detector  for  small  signals  (<-20dBm)  with 
rolloff  towards  linear  response  at  higher  input  powers.  Although  this  increases  the 
dynamic  range  of  the  measurement,  the  nonlinearity  (in  power)  must  be  corrected  for 
before  any  data  manipulation  is  done.  To  measure  this  response,  we  insert  a  known 
signal  through  a  calibrated  attenuator  and  into  the  diode.  A  generated  pulse  train  is 
then  used  to  trigger  the  signal  averager.  The  recorded  dc  waveforms  can  later  be  used 
to  determine  rf  power  input  versus  recorded  voltage.  Figure  5.4.2  shows  the  results  of 
such  a  measurement.  The  curve  in  this  figure  is  a  fifth  order  polynomial  fit  to  the  data. 
We  have  used  this  fitted  polynomial  to  smoothly  invert  recorded  voltage  waveforms  to 
rf  power  at  the  diode  input. 
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Figure  5.3.2  Amplifier 


receiver  and  data  acquisition  scheme. 
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Figure  5.4.2  RF 


5.5  Utility/Alignment  Pyroelectric  Detectors 


We  constructed  two  handheld  detectors  for  use  in  aligning  the  expanded  laser 
beam  with  the  scanner  mirrors  and  ensuring  that  the  transmitted  beam  is  striking  the 
calibration  target.  The  first  was  made  for  use  inside  the  trailer  with  an  oscilloscope  and 
is  wired  for  proportional  response.  For  light  pulses  of  duration  less  than  about  10  /xsec 
(as  in  this  case)  the  output  pulses  are  about  60  /xsec  long  with  a  peak  value  proportional 
to  the  pulse  energy  striking  the  detector.  The  detector  employed  is  a  Molectron  Pl-52 
pyroelelectric  detector  with  a  2  mm  diameter  element  and  an  Irtran  window.  The 
schematic  for  this  detector  is  shown  in  figure  5.5.1. 

The  second  detector  was  made  for  use  in  the  vicinity  of  the  calibration  target  to 
assist  in  lining  up  to  the  target.  As  oscilloscopes  are  difficult  to  handle  and  read  in  the 
field  under  daylight  conditions,  this  unit  was  made  to  be  self-contained  with  audible 
(beeper)  and  visual  (LED)  indications  of  the  incidence  of  in-band,  short  duration 
radiation  pulses.  The  schematic  for  this  unit  is  shown  in  figure  5.5.2.  The  sensor 
element  employed  is  a  Molectron  PI -71  pyroelectric  detector  with  a  1  mm  element 
diameter.  The  output  of  the  detector  is  AC  coupled  to  a  quad  comparator  to  determine 
semiquantitatively  the  pulse  energy  range.  The  comparator  outputs  are  pulse  stretched 
and  used  to  drive  a  bank  of  four  LED’s.  A  piezoelectric  buzzer  is  driven  in  parallel 
with  the  most  sensitive  channel.  External  controls  allow  for  deactivation  of  the  buzzer 
and  scaling  of  the  sensitivity  range. 
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Figure  5.5.1  Schematic  of  the  handheld  proportional  output  detector. 


Figure  5.5.2  Schematic  of  the  audio/visual  handheld  detector. 


6.  CALIBEIATION  PROCEDURES 


As  noted  in  section  3.4,  signal  returns  from  a  bl  le  polystyrene  target  were  used 
to  calibrate  the  lidar  to  a  known  reference.  We  briefly  discuss  in  this  section  the 
methods  employed  in  the  field  to  obtain  target  data  as  well  as  the  analysis  of  these 
waveforms. 

6. 1  Field  Techniques 

The  target  used  on  the  GABLE  field  measurement  was  a  six  foot  square 
mounted  at  about  a  45°  angle  to  horizontal.  (This  angle  is  adjustable  to  enable 
compensation  for  possible  elevation  differences  between  the  lidar  and  the  target.)  A 
significant  off-normal  angle  is  necessary  to  avoid  specular  reflection  effects  which  are 
exhibited  by  this  and  most  target  materials.  The  reflectivity  of  the  target  material 
measured  by  JPL  was  at  such  a  45°  angle  of  incidence. 

Initial  alignment  of  the  lidar  to  the  target  was  accomplished  by  eye  using  a 
coaligned  HeNe  laser.  Our  experience  has  shown  that  if  care  is  taken  in  this 
coalignment,  the  HeNe  and  CO2  beams  will  be  separated  by  no  more  than  about  one  to 
two  feet  over  a  path  of  up  to  3.5  km.  In  particular,  repeated  measurements  of  the 
position  of  the  HeNe  and  CO2  beams  at  the  turning  mirror  and  the  scanner  aperture 
were  found  necessary  to  ensure  a  good  alignment  after  the  telescope.  It  was  also  found 
that  verification  of  this  alignment  at  a  short  range,  about  50  to  100  meters  from  the 
trailer,  was  useful  in  reducing  the  time  needed  to  locate  the  beam  at  several  km  ranges. 
The  hand  held  pyroelectric  detector  was  used  to  locate  the  CO2  beam  at  both  the  near 
and  far  ranges. 

For  acquisition  of  the  target  data,  the  fastest  digitizing  rate  (50  nsec  per  sample) 
afforded  by  the  DSP2112FS  signal  averager  was  employed  to  obtain  sufficient  range 
resolution  to  characterize  the  return  signal  pulse.  In  addition,  we  used  an  rf  gain  about 
25  to  30  dB  lower  than  that  for  data  collection  since  the  relatively  large  signals  from 
such  hard  targets  would  otherwise  have  saturated  the  amplifier. 

Figure  6.1.1,  shows  the  recorded  signal  for  such  a  calibration  run  taken  on 
8/21/90  during  the  GABLE  program.  The  signal  amplitude  here  is  that  from  the  rf 
diode  detector.  Note  the  two  signal  artifacts  at  ostensible  ranges  of  0  and  about  0.4 
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Figure  6.1.1  Sample  target  return  data  {AU210001). 


km.  The  first  is  due  to  scatter  from  the  transmission  optics-rhomb,  telescope,  turning 
and  scanner  mirrors.  The  noise  evident  before  this  feature  results  from  RF  pickup 
from  the  laser  triggering  and  firing  circuits.  The  cause  of  the  second  artifact  is  not 
kn  jwn  but  may  be  related  to  laser  afterpulsing. 

6.2  Preliminary  Data  ReductiOii 

As  discussed  in  section  5.4,  the  recorded  signals  must  be  corrected  for  the 
nonlinearity  with  power  of  the  rf  detector.  Using  the  fit  to  the  diode  detector  response 
discussed  in  section  5.4  and  shown  in  figure  5.4.2,  the  signal  can  be  converted  from 
voltage  to  rf  power  referenced  to  the  rf  detector  input.  Figure  6.2.1  shows  such  a  plot 
of  IF  signal  power  versus  range  for  the  same  ta.get  calibration  run  as  in  figure  6.1.1. 
We  have  also  divided  out  the  measured  amplifier  gain  (preamp  and  high  gain)  so  that 
this  plot  is  of  the  IF  signal  power  at  the  output  of  the  HgCdTe  detector.  Ot^er  than  a 
change  in  scale,  the  onlv  qualitative  difference  between  the  waveforms  presented  as  rf 
power  and  those  presented  as  detector  voltage  versus  range  is  the  enhancement  in  the 
former  of  the  larger  features  due  to  compensation  for  th'.  roll-off  in  the  response  of  the 
dicde  rf  detector. 

Table  6.2.1  shows  the  results  of  this  reduction  for  12  consecutive  target 
calibration  runs  taken  on  8/21/90.  In  the  third  column,  we  have  listed  the  total  energy 
in  the  heterodyne  signal  pulse  at  the  HgCdTe  detector  output.  This  was  obtained  by 
digitally  integrating  over  the  target  return  pulse.  The  average  of  this  pulse  energy 
(referred  to  as  e^  in  sections  2  and  3)  over  these  12  runs  was  used  in  section  3.5  to 
evaluate  system  sensitivity.  The  observed  14%  va.nation  in  this  quantity  from  run  to 
run  is  greater  than  that  expected  due  to  laser  power  jitter  or  even  target  scintillation. 
(As  noted  in  section  3.5,  each  run  consists  of  the  average  of  100  laser  shots.)  It  is 
possible  that  given  the  large  range  involved  and  the  terrain  along  the  sight  path  to  the 
target  that  varying  atmospheric  conditions  such  a.s  temperature  and  humidity  were  the 
cause  of  the  observed  variation  in  the  magnitude  of  the  target  return. 

The  six  runs  listed  at  the  bottom  of  table  6.2.1  were  taken  with  a  brownish 
Styrofoam  packing  material  covering  the  target.  The  returns  from  this  material  were 
detectably  smaller  than  from  the  blue  polystyrene  foam. 
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Figure  6.2.2  shows  a  similarly  reduced  data  set  taJcen  on  8/20/90  with  the  lidar 
shooting  along  the  valley  and  out  toward  the  sea  from  Cinco  Pico  (the  GABLE  ground 
sight).  This  was  taken  concurrently  with  aircraft  measurements  along  the  valley.  The 
feature  at  about  7  km  "ange  is  a  low  lying  cloud  over  the  ocean. 


TABLE  6.2.1 


DATA  RECORD 

PEAK  POWER 
10~^  W 

INTEGRATED  RETURN 
10"^^  J 

AU210001 

7.80 

4 . 69 

AU210002 

6.66 

3.99 

AU210003 

6.81 

4 .22 

AU210004 

6.46 

3 . 87 

AU210005 

6.78 

4 . 07 

AU210006 

8.65 

5.43 

AU210007 

6.24 

3 . 78 

AU210008 

7.54 

4.48 

AU210009 

7.92 

4.77 

AU210010 

5.36 

3 . 06 

AU210011 

7.73 

4.67 

AU210012 

7.15 

4.21 

AU210013 

5.10 

3 . 00 

AU210014 

3 . 59 

2 . 03 

AU210015 

3.72 

2 . 18 

AU210016 

4.41 

2  -  58 

AU210017 

3 . 02 

1.68 

AU210018 

3 . 17 

1.78 
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